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Abstract: Wo study theoretically a multi-frequency response of electrons 
in confined graphene subject to dc-ac driven fields. We explore the pos- 
sibility for using graphene nanoribbons (GNRs) to generate and amplify 
terahertz (THz) radiations in electric field domainless regime. We discover 
two main important schemes of generation; when the frequencies are com- 
mensurate, THz generation is due to wave mixing and when they are non- 
commensurate, a single strong field suppresses space charge instability and 
any weak signals can get amplified. The use of graphene as a best substitute 
for semiconductor nanoelectronic devices is suggested. 
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1. Introduction 

Graphene is a monolayer of one atom thick with fascinating carrier trans- 
port properties. Especially, its high current density and high carrier mobil- 
ity of 44000cm^y~^s~^ [1]. But attempts to utilize these unique properties 
in graphene devices is posing some difficulties. The limitation is probably due 
to several factors including; lack of bandgap in graphene sheets, edge defects, 
disorder, among others. To overcome some of these obstacles, the dimention of 
graphene sheets has to be reduced or the geometry altered. After all, new physics 
(quantization) emerge when dimensions of materials are reduced. An infinite 2D 
graphene could become ID -I- quantization along one other direction opening 
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a gap. The resulting material is known as graphene nanoribbon (GNR). De- 
pending on the nature of the edges, one can get two symmetry groups from this 
GNR; armchair graphene nanoribbon (aGNR) or zigzag graphene nanoribbon 
(zGNR). Electron dynamics of both aGNR and zGNR have different electronic 
properties, mostly due to the berry phase and pseudospin [2]. Edge states have 
significant contribution to graphene properties, because in a nanometer size rib- 
bon, massless Dirac fermions can reach the edges within a fcmto-second before 
encountering any other lattice effects, like electron-electron interaction, electron- 
phonon interaction, etc. 

In this paper, we study the phenomenon of generating frequencies in the ter- 
ahertz (THz) range. The development of sources and sensors emitting and de- 
tecting electromagnetic waves in the terahertz regime has been the subject of 
interest for some time now. And holds great promise for graphene based THz 
metamaterials, optoelectronic devices, THz lasers, fast switching mechanisms, 
spectroscopy, wireless communication [5] . Recently, THz generations are studied 
in graphene by resonance tunneling-like configuration [6], by tunable plasmon 
excitations and light-plasmon coupling [7] and by optical pumping of graphene 
[4, 11, 8]. Bloch oscillations up to 10 THz can be generated in periodic graphene 
structures [9]. Today, semiconductor superlattices are used as sources for THz 
radiation and detection. However, GNRs are better candidates because of their 
low dimensionality, striking electronic properties and the possibility of control- 
ling these properties via applied gate voltage. Graphene is also relatively easy 
to fabricate in laboratory. 

The physical mechanism governing THz generation in graphene, when subject 
to applied electric field, can be understood in terms of ballistic trajectories of 
electrons in the quasi-momentum space. When graphene is subject to an electric 
field, ballistic acceleration of charge carriers generates to-and-fro motion of the 
whole distribution function, which varies from zero to several electron volts. It 
is a collective motion of these charges that manifest THz oscillation of carriers 
in graphene. The highly nonlinearity of graphene as a carbon allotrope and 
the fact that it has non-parabolic energy spectrum can also account for THz 
production in the material. This last effect is more applicable if the frequencies 
are commensurate. When two or more commensurate frequencies interfere in 
a region they could result in creation of fields with zero frequencies, i.e static 
fields. These bias fields are responsible for Bloch oscillations at THz frequencies 
[16]. 

The remaining of this paper is organized as follows; In section 2, we introduce the 
current density of aGNR and zGNR and imposed certain conditions to reduce 
the equations to simple forms appropriate for our systems under discussion. By 
limiting the harmonics fields to only two terms, we deduce TV characteristic 
equations for describing THz generations in section 3. The equations obtained 
in the preceding section are plotted and discussed in section 4, with conclusion 
and some recommendations for future applications in section 5. 
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2. Current density equation 

For detailed calculation of current density for GNR, see our recent paper [12]. 
To avoid book keeping, we state without repetition of the proof, the relation for 
the sheet current density when graphene is subject to an external field of the 
form E{t) = Eq + J2j EjCos(ujjt + aj) as 
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(2.1) 

Where jor is the peak current density, J„ (/?) is a Bessel function of order n and 
argument /3 = elrE/hw and = cIEq/H. For the rest of this paper we will 
consider a maximum of two harmonic frequencies Wi and 1x12 ■ Because n = 2, 
the above equation will look like 
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and gs, gv arc the spin and valley degeneracies. I ~ \/3a/2, A9 ~ its/ in + 1) for 
aGNR and I = a/2 and M = 7r(s + l/2)/(n+ 1) for zGNR. 

3. Current generation and THz amplification 

When the electric field is applied to graphene, there naturally arises two schemes 
of generation and amplification, (a) unbiased, Eq ~ Q &t even harmonics with 
commensurate frequencies and (b) biased, E ^ Q a,i non-commensurate frequen- 
cies. In the following sections, we study both scenarios in details. 

3.1. Commensurate frequencies 

If one averages out Eq.(2.2) over a period of the GNR in both sides of the 
equation, we get {jit)) = j and in the right hand side a delta function emerges 
which ensures that i>i = —^vi- Further, we consider an applied field consisting 
of purely periodic multiple harmonic frequencies, i.e uj2 = yU^i, uji = Q with 
/X being an integer or fraction. Taking the sum over 712 after linearizing with 
respect to E2 in the weak field /32 << 1 limit. This restricts the order 71,2 of 
the Bessel function to take only small values ±1. The real part of the current 
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density then becomes 



The integer /i can take even or odd values. Odd integer vahies wiU yield an imagi- 
nary current density. This means that odd harmonics are not very interesting for 
THz generations but can still exhibit NDC [13]. Substituting J(){rji2)J+i{'r(i2) ^ 
r/32 /2 in the preceding equation we arrived at 



We have obtained this equation due to wave mixing of commensurate frequen- 
cies. Direct Bloch oscillations could still be induced even for zero static field, 
but for /i equal to half-integers. We do not consider that here, so we fix /3o = 0. 
In fact, it is not hard to bring the real part of Eq.(3.2) to the form in [14], 
i.e 



However, the advantage of our equation over [14] is that /x does not have to be 
2 only. 

3.2. Non-commensurate frequencies 

If the frequencies are not commensurate, then the condition vi + (11^2 = is lifted 
and llI2 7^ i'lji. However, it is still possible to amplify frequencies in the THz 
domain if one considers wi and uj2 as pump and probe frequencies respectively. 
Both frequencies belong to THz range, i.e ujit > 1 and UJ2T > 1. This approach 
has been adopted in [15, 16] for superlattices and in experiment for generating 
THz using resonance tunneling-like configuration in graphene [6] and amplifying 
small frequencies in epitaxially grown graphene heterostructures [8]. Here, the 
THz is generated because the pump field excites electrons, the space charge 
instability is suppressed by the strong (pump) field while the week probe signals 
get amplified. 

The real part of Eq.(2.2) for the two frequencies which are not related (ai_2 = 0) 
takes the form 
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where we have assume that vi = 1^2 = 0. 

Most often, we shall be using the ratio £1^2/ Ecr which is defined as 

-fT^ = /3i.2Wi,2T, with Ecr = -7- (3.5) 

her S.Lt 

instead of just /3i,2- 

4. Results and Discussion 

In Fig.l, we have demonstrated dependence of non-linear current density on 
the harmonic index ^ for week (top) and strong (buttom) ac amplitudes. In 
both cases, some ^ values give positive (j_) and negative (j+) current density. 
The series resulting in j+ yields Bloch oscillations of the current density that 
decays faster to j = from above as in Fig. 2 (right). The series that results 
in j_ produces Bloch oscillations that is asymptote to j = from below as in 
Fig. 2 (left). This means electronic oscillations persist for some time before dying 
off and thus the ^-series that give j_ is the better option for THz production. 
They arc two sub-categories of j±. j+: ^lodd = 1,9,..., ^JLeven — 4,6,8,10... 
at low ac amplitudes and ^odd = 1,5,..., fieven = 4, 8, . . . at high amplitudes. 
j_: fiodd = 3, 5, 7, . . ., fieven = 2, ... at low ac amplitude and ^J,odd = 3, 7, 9, . . ., 
Meuen = 2, 6, 10, . . . at high amplitudes. From the graphs, we found that fieven = 
2,4 and fiodd = 1,3 are robust against the ac field amplitudes, Ei, E2. 




Fig 1. Non-linear current density with harmonic index, fi. cos{a) = 1 
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The graphs in Fig. 2 show the behavior of terahertz eurrent density on the ac 
field amphtudes for armchair and zigzag graphene nanoribbons at /i = 2,4. 
For low Ei/Ecr values at fir = 2, the absolute current steadily increases to 
maximum before falling. The curve then begins to oscillate above certain ac 
threshold amplitude, Ei^nin- Under our default parameters, Ei,„iin = 2iJcr for 
aGNR and Ei_„iin = 2\f^Ecr for zGNR. This oscillations are predicted to lie 
within THz range [14]. 
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Fig 2, THz current density against ac field amplitudes for armchair and zigzag graphene 
nanoribbons. (Left) second harmonic, fi = 2 and (right) fourth harmonic fi = 4. cos{a) = 1. 



We also studied in Fig. 3 the combined effect of Thz current with reduced ac 
amplitudes and frequencies. Oscillations disappear when fir << 1. This region 
is not feasible for THZ generation. However, for fir > 1 oscillaions are more 
pronounced and a graphene device under this condition can operate effectively 
to produce THz frequencies. 

Within the scheme of commensur ability of frequencies, we have finally demon- 
strate the behavior of THz current on both phase difference and ac amplitudes 
in Fig. 4. There are some phase differences where there is no oscillations, this is 
indicated on the contour as straight lines. Interestingly, current peaks when the 
two ac fields are out of phase, i.e a = tt and Ei = 5Ecr- 

Now, we turn to the case of non-commensurate frequencies. In Fig. 5, we plot 
THz current with ac amplitudes in the presence of another stronger ac amplitude 
but with weak frequency and biased static field. It is this weak frequency that 
excites electrons after the electric charge instability is suppressed by the strong 
ac field amplitude. The excited electrons Bloch oscillates at rather magnified 
frequency within the THz range. However, the nature of oscillations in the figure 
is different from what we have seen earlier. Here, the curve falls quickly and then 
begins to oscillate almost immediately. 
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Fig 3. THz current density against ac field amplitudes for armchair and zigzag nanoribbon 
at various Or values, cos (a) = 1 
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Fig 4. Three dimensional plot showing THz current oscillations at varying ac field and phase 
difference. 



The graph in Fig. 6 also shows effects of both the ac fields Ei, E2 on THz current. 



Rabiu Musah et al. /Terahertz generation and amplification in graphene nanoribbons 8 




0.1 i 




Fig 6. Three dimensional plot showing THz current oscillations at varying ac fields. L02 = 
0.7a;i, Eo = 2.5Ecr 



5. Conclusion 



We have studied theoretically, two schemes of THz production; generation and 
amplification at commensurate and non-comnicnsurate frequencies respectively. 
For commensurable frequencies, we discovered two harmonic series, fJ-odd = 
2,4,... and fj-even = 1,3,... which are robust against high ac field amplitudes 
and at which values production is suitable. Generation at the commensurate 
frequencies is due to wave mixing at zero bias field. For incommensurable fre- 
quencies, there is THz amplification of small ac signals. A strong ac amplitude 
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is ncccessary to suppress charge instability at non-zero bias field. 
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